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Abstract 

A hybrid Finite Element Method (FEM)/Method of Moments (MoM) technique in 
conjunction with the Asymptotic Waveform Evaluation (AWE) technique is applied to obtain 
radar cross section (RCS) of a cavity-backed aperture in an infinite ground plane over a 
frequency range. The hybrid FEM/MoM technique when applied to the cavity-backed aperture 
results in an integro-differential equation with electric field as the unknown variable, the electric 
field obtained from the solution of the integro-differential equation is expanded in Taylor series. 
The coefficients of the Taylor series are obtained using the frequency derivatives of the integro- 
differential equation formed by the hybrid FEM/MoM technique. The series is then matched via 

the Pade approximation to a rational polynomial, which can be used to extrapolate the electric 
field over a frequency range. The RCS of the cavity-backed aperture is calculated using the 
electric field at different frequencies. Numerical results for a rectangular cavity, a circular 
cavity, and a material filled cavity are presented over a frequency range. Good agreement 
between AWE and the exact solution over the frequency range is obtained. 
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1. Introduction 


The electromagnetic characterization of cavity backed apertures is of importance in 
understanding the scattering properties and in electromagnetic penetration/coupling studies. 
Recently, there has been a considerable interest in analyzing cavity backed apertures in an 
infinite ground plane. Various analytical and numerical techniques have been applied for two 
dimensional cavity backed apertures [1-6]. For three dimensional problems, mode matching has 
been used for rectangular [7] and spherical [8] cavity backed apertures. A Method of Moments 
(MoM)/modal [9] approach is used recently to analyze apertures formed by a rectangular cavity 
recessed in a ground plane. These methods are restricted to cavities with regular shapes, where 
fields can be written in modal form. In [10], a boundary integral method is used to analyze the 
scattering from three dimensional cavities via a connection scheme. Though useful for savings 
in computer memory, this method leads to dense matrices. Also, the accumulation errors due to 
the connection algorithm are not negligible as the number of subsections increase. In the case of 
deep cavities, high frequency techniques such as those proposed in [11] and [12] could be 
effectively implemented. Unfortunately, these techniques are not suitable when the cavity is 
filled with inhomogeneous materials. In [13], a Finite Difference Time Domain (FDTD) method 
is applied for large structures. However, the method of FDTD sometimes results in inaccurate 
results due to differencing, staircasing and dispersion. Jin and Volakis [14] used a finite element- 
boundary integral formulation which employs the boundary integral equation (or Method of 
Moments-MoM) to formulate the fields external to the cavity accurately. 

In the combined FEM/MoM technique, FEM is used in the cavity volume to compute the 
electric field, whereas MoM is used to compute the magnetic current at the aperture. For the 
combined FEM/MoM technique, the cavity is divided into tetrahedral elements and the aperture 
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is discretized by triangles. Using the Galerkin’s technique, simultaneous equations are generated 
over the subdomains and are added to form a global matrix equation. This results in a partly 
sparse and partly dense, symmetric complex matrix, which can be solved either by a direct 
solver or by an iterative solver. The electric field hence obtained is used to compute the RCS of 
the cavity-backed aperture. 

To obtain RCS over a range of frequencies using the combined FEM/MoM technique, one 
has to repeat the calculations over the frequency range of interest. If the RCS is highly frequency 
dependent, one needs to do the calculations at the finer increments of frequency to get the 
accurate representation of the frequency response. This can be computationally intensive and for 
an electrically large cavity with electrically large aperture, it can be computationally prohibitive 
despite the increased power of the present generation of computers. To alleviate the above 
problems, the application of Asymptotic Waveform Evaluation (AWE) has been proposed for the 
timing analysis of VLSI circuits [15]. The AWE technique is finding increasing interest in 
electromagnetic analysis of microwave circuits [16-18]. Recently a detailed description of AWE 
applied to frequency domain electromagnetic analysis is presented in [19]. AWE has been used to 
predict RCS of PEC bodies over a frequency range [20] and input characteristics of a 
cavity-backed aperture antenna over a frequency range [21]. 

In this report, the application of AWE for predicting the RCS over a range of frequencies 
for a cavity-backed aperture using a combined FEM/MoM technique is described. In the AWE 
technique, the electric field is expanded in a Taylor series around a frequency. The coefficients of 
the Taylor series (called ‘moments’) are evaluated using the frequency derivatives of the 
combined FEM/MoM equation. These moments are then matched via the Pade approximation to 
a rational polynomial. Using the rational polynomial, the electric field distribution in the cavity 
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can be obtained over a frequency range. Using this field distribution, the RCS of the 
cavity-backed aperture in an infinte ground plane is calculated at any frequency within the 
frequency range. 

The rest of the report is organized as described below. A brief description of hybrid 
FEM/MoM technique to calculate RCS from a cavity-backed aperture is given in section 2. In 
section 3, AWE implementation for the combined FEM/MoM technique is described. Numerical 
results for a rectangular cavity, a circular cavity, and a material filled cavity are presented in 
section 4. The numerical data are compared with the exact solution (calculated at each frequency 
using the hybrid FEM/MoM technique) over the frequency range. CPU time and storage 
requirements for AWE formulation are given for each example and are compared with those 
required for exact solution at each frequency. Concluding remarks on the AWE technique are 
presented in section 5. 
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2. Hybrid FEM/MoM Technique 

The geometry of the problem to be analyzed is shown in figure 1. For linear, isotropic, and 
source free region, the electric field satisfies the vector wave equation: 

Vx^VxeJ -& 2 £ r E = 0 (1) 

where p r , E r are the relative permeability and relative permittivity of the medium in the cavity 
and k is the free space wavenumber. The time variation exp (/tot) is assumed and suppressed 
throughout this paper. The electric field E is solved via hybrid FEM/MoM technique assuming 
that cavity-backed aperture is illuminated by a harmonic plane wave, H jnc . Following the proce- 
dure described in [14], equation (1) can be written as 

JJJ (VxT) .(’IvxEW^EjJjT.Erfv-ymnJ j (Txft) 

V vH r y v s ap 

= J (T x n) • H inc ds (2) 

S a, 

where T is the vector testing function. S ap is the aperture surface (see figure 1). is the 

scattered magnetic field and n = z at S . 

The volume of the cavity is subdivided into small volume tetrahedral elements. The 
electric field is expressed in terms of the edge vector basis functions [22], which enforce the 
divergenceless condition of the electric field explicitly. The vector testing function is also 
expressed in terms of the edge vector basis functions following the Galerkin’s method. The 
discretization of the cavity volume into tetrahedral elements automatically results in discretization 
of the surface S ap into triangular elements. The volume and surface integrals in equation (2) are 
carried out over each element to form element matrices, which are assembled to form global 
matrices. Equation (2) can be written in matrix form as 
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A (it) e(k) = b (it) 


( 3 ) 


A (it) is a partly sparse, partly dense complex symmetric matrix, b(k) is the excitation vector, and 
e(k) is the unknown electric field coefficient vector. A(k) is evaluated as a sum of three matrices. 


where 


A (it) = Aj(it) +A 2 (k) +A 3 (k) 


A,(*) = JJJ (VxT) .flvxEWk 2 e r JJjT.Edv 

V V v 




S V s 

ap \ ap 


ds 


A 3 ( k) = ijj (V.T) (V'.M)^UhEl ds '\ ds 


R 


b(k) = 2wJJ (Txn) *H inc ds 


( 4 ) 

( 5 ) 

(6) 

( 7 ) 

( 8 ) 


Equations (6) and (7) are obtained by making use of the equivalence principle and image theory 
[23] and follow the procedure given in [24]. T 5 = Txfi and R is the distance between source 

point and the observation point. M is the equivalent magnetic current over the aperture S ap . V' 

indicates del operation over the source coordinates and ds' indicates the surface integration over 
the source region. Equation (8) is calculated assuming a harmonic plane wave 

H ( „ c = (* H xl + 9H yl + tHJe~ ,% ” = + (9) 

E iTC = 1 „H,.„ c xk,. (10) 


where 
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= -k o [ksinG^cosfy + ^sinG.sin^. + fccosGJ 

(11) 

H xi - ( sin a cos 0. cos 4^ + cosasin^.) /r\ o 

(12) 

H yi = (sinacosGj.suKjK- cos a cos fy) /'ll o 

(13) 

H zi = (-sinasinG^ /r\ a 

(14) 

H 0i = | H ,*c| sina 

(15) 

= | H «»c| coscx 

(16) 


in which T] is the free space intrinsic impedance and a represents the polarization angle of the 
incident field. When a = 0, then H zi = 0 which corresponds to H-polarization and when 
a = Jt/2 , then E zi = 0 which corresponds to E-polarization. 

The matrix equation (3) is solved at any specific frequency, f g (with wavenumber k g ) 
either by a direct method or by an iterative method. The solution of the equation (3) gives the 
unknown electric field coefficients which are used to obtain the electric field distribution. Once 
the electric field E is found and hence the magnetic current M on the aperture, the far zone 
scattered field can be computed. 


H - <r) L. - -f ifJI (§§+W) ,MU - 3 ' )e 

*o s 


7'Jt o sin0 (atcosi}) + vsirx))) 


dxdy 


(17) 


where ( r, 0, <j)) are the usual spherical coordinates of the observation point. The scattering cross 
section is then given by 


c 


lim 4rcr 2 

r — » oo 


|H, e „,(r)l 2 

|H,„ c <r)f 


(18) 
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where 


(19) 


H ,„,( r )| 2 = |// e ,| 2 + KJ 2 

| H ,„( r )| 2 = | H ei | 2 + |^,| 2 ( 20 ) 


3. AWE Implementation 

The RCS given in equation (18) is calculated at one frequency. If one needs RCS over a 
frequency range, this calculation is to be repeated at different frequency values. Instead, AWE can 
be applied to obtain the frequency response over a frequency range. The general implementation 
of AWE for any frequency domain technique used for electromagnetic analysis is given in detail 
in [19]. The solution of equation (3) gives the unknown electric field coefficient vector e (k o ) at a 
particular frequency / . Instead e (k) can be expanded in Taylor series as 


e(k) = £ m n (k-k o ) n 


( 21 ) 


71 =0 


with the moments m n given by [19] 


m n = A ' (*„) 


b in) (k c ) " 

2-t n \ 


n! 


( 22 ) 


<7 = 0 


A <<?) (k o ) is the derivative with respect to k of A(k) given in equation (4) and evaluated at k g . 
Similarly, b (k o ) is the <yth derivative with respect to k of b(k) given in equation (9) and 
evaluated at k . The Kronecker delta is defined as 

o go 



i 

o 


q = 0 

q h 0 


( 23 ) 


The <7* derivatives of A(k) and b(k) are evaluated and are given in detail in the Appendix. 
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In many cases, the Taylor series expansion gives fairly good results. However, the 
accuracy of the Taylor series is limited by the radius of convergence. It will not converge to the 
right answer beyond the radius of convergence, and it sometimes requires a large number of terms 
to converge over a frequency range. In such cases, one may want to replace Taylor series 
expansion with a rational function called Pade approximation [15] to improve the accuracy of 
the numerical solution. 

To obtain Pade approximation, the Taylor series expansion in equation (21) is matched 
with a rational polynomial [15] 


X m n( k - k o ) 

rc = 0 




(24) 


where 

P L (^ — = — kp) ^2 ( k — k o ) & L ( k ~ k o ^ 

and 

7 M 

8* (*-*„) = b 0 + b,(k-kj + b 2 (k-kj + + b u (k-k 0 ) 

b o is set to 1 as the rational function can be divided by an arbitrary constant. Since there are 
(L+M+l) unknowns, (L+M) moments of the Taylor series should be matched. Equating the 

coefficients for powers (k-k 0 ) L+i (k-k o ) L + M , the coefficients of Q M (k- k g ) can be 

obtained solving the matrix equation 


m L-M + 1 m L-M + 2 m L 

m L-M + 2 m L-M+ 3 •” m L + 1 


' b M 

b M-\ 


m L + 1 
m L + 2 

m L m L+ 1 m L + M- 1 . 


. b 1 . 


™L + M_ 
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The numerator coefficients can be found by equating the powers (k-k 0 )° (k - k o ) L 

a o = m o 

a \ = m l +b l m o 

a 2 ~ m 2 + b\ m \ + b 2 m 0 

min (L, M) 

a L = m L+ Z b i m L-i 
i = 1 

Once the coefficients of the rational polynomial are obtained, equation (21) can be 
rewritten as 

a + a , (k-k ) +a 7 ( k-k ) 2 + +a, (k-k ) L 

e (k) = -2 — — ^ 4 i- 24 (26) 

1+*, (*-*„) +b 2 (k-k a ) + +b M (k-k 0 )" 

For a given amount of computational effort, one can easily construct a rational 
approximation which has smaller error than a polynomial approximation. Also for a fixed value of 
L+M, the error is smallest when L-M or L=M+7[17], Using equation (26), the electric field 
coefficients at frequencies around the expansion frequency are obtained. The electric field hence 
obtained is used to compute the scattered magnetic field given in equation (17) and finally the 
backscattering cross section using equation (18). 

4. Numerical Results 

To validate the analysis presented in the previous sections, a few examples are considered. 
RCS calculations over a frequency range are done for a rectangular cavity, a circular cavity, and a 
cavity filled with lossy material. The numerical data obtained using AWE are compared with the 
results calculated at each frequency using the computer code CBS3DS [25], which implements 
the combined FEM/MoM technique [14]. We will refer to the latter method as “exact solution.” 
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From section 3, it can be observed that the inverse of matrix A (k o ) is found once and is used 
repeatedly to find AWE moments. Due to the hybrid FEM/MoM technique, matrix A (k o ) is 
partly sparse and partly dense. The Complex Vector Sparse Solver (CVSS) [26] is used to LU 
factor the matrix A (i t ) once and the moments are obtained by backsolving equation (21) with 
multiple righthand sides. All the computations reported below are done on a SGI -Indigo2 (with 
150MHz, IP22 processor) computer. 

(a) Square Cavity : A square cavity in an infinite ground plane is considered (fig. 2 with 
a=lcm, b-lcm, and c=2cm). Backscattering calculations are done with an incident angle 

0 = 0° and <() = 0° . The discretization of the square cavity resulted in 3590 total unknowns, 
and the order of the dense matrix due to MoM is 133. Figure 3a shows the radar cross section 
over the frequency range 15GHz to 25GHz, calculated using Taylor series for a H-polarized 
incident wave. The Taylor series moments are calculated at 20GHz. Figure 3b shows the radar 

cross section over the frequency range 15GHz to 25 GHz calculated using Pade approximation. 
It can be seen from Figure 2b that Taylor series gave good results over the frequency range 
18GHz to 22GHz. Beyond this frequency range, there is no improvement in accuracy, even by 

adding more terms to the Taylor series. However, figure 3b indicates that Pade approximation 
gave good results over the frequency range 15GHz to 25GHz with L=5 and M=5, and well 
behaved convergence is observed with increase in the orders of numerator and denominator of 

Pade approximation. The timings for the calculations performed using CBS3DS and Pade 
approximations are given in Table 1. Note that the tunings for Taylor series expansion and Pade 
approximation are the same except for the cost of calculating Pade coefficients from the Taylor 
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series coefficients. It is observed that the cost of generating moments and the Pade coefficients 
is very minute compared to the cost of matrix generation and solution. 

(b) Circular Cavity: A circular cavity (radius -0.305cm and height=0.3cm) in an infinite ground 
plane is considered (fig. 4a). The discretization of the circular cavity resulted in 1327 total 
unknowns, and the order of the dense matrix due to MoM is 132. Backscattering from this 

cavity is calculated over the frequency range 10GHz to 50GHz with the incident angle 0 = 0 

and (j> = 0° and plotted in Figure 4b. The Taylor series moments are calculated at 30GHz. 

Taylor series expansion is calculated with five moments and the Pade approximation is 

calculated with L=3 and M-2. It can be seen that Pade approximation could give accurate 
results throughout the frequency range, whereas Taylor series gave accurate results only within 
the frequency range 21 GHz to 39GHz. The timings for calculation of backscattering cross 

section using “exact method” and the Pade approximation are given in Table 1. 

(c) Rectangular cavity with lossy material: A rectangualr cavity is considered as another 
example (fig. 2 with a=lcm, b=0.25cm, and c=0.25cm). The cavity is filled with lossy material 
with dielectric constants e r = 2.2 -y'1.5 and |l r = 1.8 -jOA. The cavity is discretized using 

tetrahedral elements resulting in 3218 unknowns. The order of the dense matrix due to MoM is 
275. The backscatter cross section is calculated over the frequency range 10GHz to 50GHz with 

the incident angle 0 = 0° and <|> = 0° and plotted in Figure 5a. The incident wave is assumed 
to be H-polarized. Taylor series expansion is calculated with five moments at 30GHz and 

accordingly the Pade approximation is calculated with L=3 and M=2. It can be seen that Pade 
approximation gave accurate results throughout the frequency range, whereas Taylor series gave 
accurate results only within the frequency range 22GHz and 38GHz. The backscattering cross 
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section calculations are also carried out for an E-polarized incident wave with incident angle 

0 = 0° and <j> = 0° and plotted in Figure 5b. Pade approximation is calculated with L=5 and 
M-4. Taylor series approximation is calculated with nine moments at 30GHz. It can be seen 

even for this case that Pade approximation gave accurate results throughout the frequency range 
10GHz to 50GHz, whereas Taylor series is accurate only within the frequency range 20GHz to 
38GHz. The timings for calculation of backscattering cross section using “exact method” and 

the Pade approximation are given in Table 1. The discretization for the “exact” calculations 
using CBS3DS at frequency points beyond 40GHz resulted in 5848 unknowns and order of the 

dense matrix due to MoM is 421. It can be noted that both Taylor series and Pade 
approximation are evaluated at 30GHz and hence the dicretization that is used at 30GHz is 
accurate enough to calculate the hybrid FEM-MoM matrix and the derivative matrices. The 

results obtained using Pade approximation show a good agreement with the “exact” calculations 
with denser gridding even beyond 40GHz. 

Comment on Storage : In all the above examples, when solving a matrix equation, one needs to 
store the matrix A (k 0 ) for exact solution at each frequency. For « th order AWE, one needs to 

store n number of matrices (k o ) , q=l,2,3,...n), along with the matrix A (k o ) . For 

electrically large problems, this could impose a burden on computer resources. This problem can 
be overcome by storing the derivative matrices, A ^ (k o ) out-of-core, as the derivative matrices 
are required only for matrix-vector multiplication. 
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Table 1: Comparison of CPU timings for the numerical examples presented in Sections 4a, b 
and c 


Problem 

Method 

Matrix Fill 

(secs) 

LU Factor 

(secs) 

Total 

Time 

(secs) 

(a) Rectangular 
cavity 

(a=Jcm, b -0.5 cm, 
c=0.5cm) 

CBS3DS 
(13 freq. points) 

2081.7 

312.0 

2393.7 

Pad e'(L=5, 
M=5) 

(100 freq. points) 

354.7 

25.0 

379.7 

(b) Circular cavity 
( radius = 0. 305cm, 
height=0.3cm) 

CBS3DS 
(41 freq. points) 

6342.7 

270.60 

6613.3 

Pade(JL=5,Af=2) 
(400 freq. points) 

323.7 

6.64 

330.34 

(c) Material filled 
rectangular 
Cavity 

(a=lcm,b=0.25cm, 

c=0.25cm) 

CBS3DS 
(41 freq. points) 

19526.9 

( 1 0GHz~40GHz) 
+ 

14880.0 

(41GHz-50GHz) 

1110.73 

(10GHz-40GHz) 

+ 

1413.5 

(41GHz-50GHz) 

37021.13 

Pade(L=5,A/=2) 
(400 freq. points) 

(H-Pol) 

1330.1 

34.67 

1364.77 

Pade'(I=J,M=4) 
(400 freq. points) 

(E-Pol) 

1429.21 

34.67 

1463.88 
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5. Concluding Remarks 


The AWE technique is applied to the hybrid FEM/MoM technique to calculate the radar 
cross section of a cavity-backed aperture over a frequency range. Examples of a rectangular 
cavity, a circular cavity, and a material filled cavity are considered to validate the analysis. Both 

Taylor series approximation and Pade approximation are calculated for all the examples. It can 

be noted that for the same computational effort, the Pade approximation proved to be superior 
in terms of wider bandwidth. Timing comparisons are done for calculating radar cross section 
over a frequency range using AWE and using ‘exact’ calculation at each frequency point. AWE 
is found to be superior in terms of the CPU time. It may be noted that although calculations are 
done in frequency increments of 0.1 GHz for eamples presented in this paper, frequency 
response at even finer increments can also be calculated at a very nominal cost. This is 
particularly important when there are sharp nulls present in the frequency response. 
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Appendix 

Derivatives oiA(k) and b(k) w.r.t. k 


The frequency derivatives of A(k) and b(k) are evaluated and are given below. From 
equation (4): 


A {q) (it) = — (k) = A\ q) (it) +A ( 2 q) ( k ) +A (q) ( k ) +A iq) ( k ) q=0,l,2,3, (A.l) 

dk q 

From equation (5) 

A, (0) (it) = JJJ l(VxT) • (VxE)dv-it 2 e r JJjT*EJv (A.2) 

V V 

A, (1) (it) = -2iteJJjT«EJv (A.3) 

v 

a[ 2) (it) = -2e r JJ|T • Edv (A.4) 

v 

Aj ?) (k) = 0 qZ3 (A.5) 



A ( 2 q) <*) 


j j T s • ( j j m( 3- ) (-y/?) «■ " 3 + 2qk HR ) 9 " ' 2 + k 2 (,-jR) 9 ' ’] exp (-jkR) ds' ds far q> 1 (A. 8) 

S afl ( S ap J 
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From equation (7) 


(*> = j (V«r.) |j f (V’ • M) ZZU^ dsjds 

*"* ap ^ ap 

A™ ( k ) = J | (V«T 5 ) { J J (V' • M) (-./*) {q ~ ° exp (-./**) *'} * 


From equation (8) 


b (0) ( k ) = 2 /topJ J (Txn) •H |nc ds 


Noting from equation (9) that 


H. 


jk(x t +y x + z,) 

He 


where 


H. = ZH xi + fH yi + zH zi 
= jcsin0.cos<j) ( . 

>>! = }>sin0.sin<t> ( 

and z, = zcosSj. 


Equation (A.l 1) can be rewritten as 




jk(x , +y t + z,) 


(A.9) 
(A. 10) 

(A. 11) 

(A. 12) 


(A. 13) 
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For < 7^0 


b w (k) = ‘ f ' b(k) - 

dk q 


= 2# 0) ?T I 0 J J (Txfi) • +^j +Zj) 


( 9 _I) jr* (jc t +>-j + z,) 


ds 


+ 


S„, 


(Txn) • H f (x, +>, +Zj) 


<7 


jk (x, +>, +Z!) 

e 


ds 


(A. 14) 
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Figure 1 Geometry of a cavity-backed arbitrarily shaped aperture in an infinite ground plane. 
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Fig. 2 Geometry of the rectangular cavity backed aperture. 
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(dB) 
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Figure 3a Frequency response calculation for the rectangular cavity shown in figure 2 (a=lcm, 
b=lcm, c=2cm, e r =1.0, \l r =1.0) using Taylor series approximation. 
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Frequency(GHz) 


Figure 3b Frequency response calculation for the rectangular cavity shown in figure 2 (a=lcm, 
b=lcm, c=2cm, e r =1.0, \L r =1.0) using Pade approximation. 
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Figure 4. Frequency response calculations of back scattering from an air-filled circular cavity, 
(a) Geometry of the circular cavity (b) Backscattering cross section versus frequency. 
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Figure 5a Frequency response calculation for the rectangular cavity shown in figure 2 ( a=lcm, 
b=0.25cm, c=0.25cm, e r =2.2-jl.5, \i r =1.8-jO.I). H-Polarization 
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Figure 5b Frequency response calculation for the rectangular cavity shown in figure 2 (a=lcm, 
b=0.25cm, c=0.25cm, z r =2.2-jL5, \i r =1.8-j0.1). E-Polarization. 
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